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I. INTRODUCTION 

During the last decade, regular arrays of magnetic 
nanoparticles have been deeply investigated. Besides the 
basic scientific interest in the magnetic properties of these 
systems, there is evidence that they might be used in 
the production of new magnetic devices, 'l', 2] Different 
geometries have been considered, including dots, rings, 
tubes and wires. Recent studies on such structures have 
been carried out with the aim of determining the stable 
magnetized state as a function of the geometry of the 
particles. 0, 0, B 0| In particular, the study of highly 
ordered arrays of magnetic wires with diameters typically 
in the range of tens to hundred nanometers is a topic of 
growing interest. 0, d, S This is a consequence of 
the development of experimental techniques that lead to 
fabricate in a controllable and ordered way such arrays. 
(Til [T^ The high ordering in the array, together with the 
magnetic nature of nanowires, give rise to outstanding 
cooperative properties of fundamental and technological 
interest. [31 

Bistable nanowires are characterized by squared- 
shaped hysteresis loops defined by the abrupt reversal of 
the magnetization between two stable remanent states. 
[13, [iHl In such systems, effects of interparticle interac- 
tions are in general complicated by the fact that the dipo- 
lar fields depend upon the magnetization state of each 
element, which in turn depends upon the fields due to 
adjacent elements. Therefore, the modelling of interact- 
ing arrays of nanowires is often subject to strong sim- 
plifications like, for example, modelling the wire using a 
one-dimensional modified classical Ising model. 0, [3| 
Zhan et al. [l9| used the dipole approximation includ- 



ing additionally a length correction. J. Velazquez and M. 
Vazquez ,17, 18,J considered each microwire as a dipole, in 
a way that the axial field generated by a microwire is pro- 
portional to the magnetization of the microwire. Never- 
theless, this model is merely phenomenological since the 
comparison of experimental results with a strictly dipo- 
lar model shows that the interaction in the actual case is 
more intense. They have also calculated the dipolar field 
created by a cylinder and expanded the field in multipo- 
lar terms, (20| showing that the non-dipolar contributions 
of the field are non negligible for distances considered in 
experiments. In spite of the extended use of the dipole- 
dipole approximation, a detailed calculation of the valid- 
ity of approximations on the dipolar energy has not been 
presented yet. Also micromagnetic calculations [2l|, [S^l 
and Monte Carlo simulations [1^ have been developed. 
However, these two methods permit to investigate arrays 
with just a few wires. In this paper we investigate the 
validity of the dipole-dipole approximation and include 
additional terms that lead us to consider large arrays and 
the shape anisotropy of each wire. 



II. CONTINUOUS MAGNETIZATION MODEL 

Geometrically, nanowires are characterized by their ra- 
dius, R, and length, L. The description of an array of 
N wires based on the investigation of the behavior of 
individual magnetic moments becomes numerically pro- 
hibitive. In order to circumvent this problem we use a 
continuous approach and adopt a simplified description 
in which the discrete distribution of magnetic moments 
in each wire is replaced with a continuous one, defined 



2 



by a function M(r) such that 'M.{r)5V gives the total 
magnetic moment within the element of volume SV cen- 
tered at r. We recall that Etot is generally given by the 
sum of three terms corresponding to the magnetostatic, 
Edip, the exchange, E^x, and the anisotropy contribu- 
tions. Here we are interested in soft or polycrystalline 
magnetic materials, in which case the anisotropy is usu- 
ally disregarded, [g] 

The total magnetization can be written as M(r) — 
^^-^Mi(r), where Mi(r) is the magnetization of the 
i-th nanowire. In this case, the magnetostatic poten- 
tial U{r) splits up into N components, Ui{r), associ- 
ated with the magnetization of each individual nanowire. 
Then, the total dipolar energy can be written as Edip = 

J2f=i Edip{i) + Y.^Ji Eintihj), where 



(1) 



1. Self energy of a nanowire 

The reduced dipolar self-energy has been calculated by 
Landeros et al. 5] and takes the form 



(4) 



where F2i[a:] — F2i[— 1/2, 1/2, 2, a;] is a hypergeometric 
function. Note that in Eq. Q the energy of each wire 
depends only on the ratio R/ L. 



2. Interwire magnetostatic coupling 

The interaction between two nanowires is obtained us- 
ing the magnetostatic field experienced by one of the 
wires due to the other. Details of these calculations are 
included in Appendix A, giving 



is the dipolar contribution to the self-energy of nanowire 
i-th, and 



Ernt = Mo / M,(r)VC/j (r) dv 



is the dipolar interaction between them. In the dipolar 
contribution to the self-energy an additive term indepen- 
dent of the configuration has been left out. [l^l 

In this work we investigate bi-stable nanowires in 
which case, E^x = Yli'i=iEex{i) — 0. On the basis 
of this result, the total energy of the array can be writ- 
ten as 



N 



N-1 N 



Et. 



Y,Eself{i)+Y. E (3) 



i=l 



i—1 



where Eseif{i) — Edip{i) is the dipolar self-energy of each 
wire, and Eint(i,j) is the dipolar interaction energy be- 
tween wires i-th and j-th. 



A. Total energy calculation 

We now proceed to the calculation of the energy terms 
in the expression for Etot- Results will be given in units 
of HoM^V, i.e. E = E/^iqM^V, where V = ttR^L is 
the volume of the nanowire and Mq is the saturation 
magnetization. 

In order to evaluate the total energy, it is necessary to 
specify the functional form of the magnetization for each 
nanowire. We consider wires with an axial magnetization 
defined by Mi(r) = MoCTiZ, where z is the unit vector 
parallel to the axis of the nanowire and Gi takes the values 
±1, allowing the wire i to point up ((7^ = +1) or down 
(tTi = —1) along z. 



£;.„,(^, j) = 2a.a, / ^ Jo Jf i^f] (1 - e"^) 



(5) 

(2) where Jp is a Bessel function of first kind and p or- 
der and Sij is the center-to-center distance between the 
magnetic nanowires i and j. The previous equation al- 
lows us to write the interaction energy of two wires as 
Eint[i,j) = aia-jEmtiSij) = ±EtntiSij), where the sign 
+ (— ) corresponds to ai — aj [cTi ^ aj), respectively. 



B. Results 

1. Two wires system 

The general expression giving the interaction between 
wires with axial magnetization is giving by Eq. ([5]). This 
expression has to be solved numerically. However, wires 
that motivate this work 

0, i, i, satisfy L/R > 1, 
leading us to expand Ji as 



Ji{x) 



X X 



2 16 ^ kW{k + 2) 



Then we can approximate Eq. (O by 

?2 °° 



Eint{Sij) 
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(7) 



where A indicates the order of the expansion. As an il- 
lustration, the first and second terms in the sum are 



9i 
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1 

ai ■ 



(8) 



and 
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R^ 



ai 



^d^^] (9) 
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where ai = .^l + LV^^-, ^2 = (1 - 2L^ISl)/a\ and 

as = (1 + SiVSSi'j - 8L^/Sfj)/al. Figure 1 illustrates 
the interaction energy between two identical nanowires 
with parallel axial magnetization as a function of 2R/Sij. 
When the two wires are in contact, 2R/Sij = 1; when 
they are infinitely separated, 2R/Sij = 0. In this figure 
the solid line represents the numerical integration of the 
interaction energy, Eq. ([5]), the dashed line is given by 
the first order approximation of this energy, Eq. ([7]) and 
the dotted line corresponds to the second order approx- 
imation. From this figure we observe that a first order 
approximation gives a reasonable approach to Eq. ([5]) 
for L/R ^ 1. We can conclude that the first term in the 
expansion in Eq. ([7]) gives a very good approach to Eq. 
© for 2R/S^J < 0.6, and L/R > 1. 




2R/S.. 
ij 



the magnetostatic interaction, Eq. (O, and the dipole- 
dipole approximation, Eq. (jlOp . between two identical 
nanowires as a function of 2R/Sij. These results are 
illustrated in Figure 2 and lead us to conclude that the 
dipole-dipole approximation overestimates the real inter- 
action, except for very apart wires, which is not the usual 
case in an array. Also for large aspect ratio wires the ap- 
proximation becomes worst. 
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FIG. 2: Ratio between the magnetostatic interaction, Eq. (5), 
and the dipole-dipole interaction, Eq. (10), as a function of 
2R/Sij. Different aspect ratio L/R are considered. 



FIG. 1: Interaction energy between two identical nanowires 
with parallel axial magnetization. The solid line corresponds 
to the numerical integration of Eq. (5), the dashed line cor- 
responds to the first order approximation of Eq. (7) and the 
dotted line corresponds to the second order approximation of 
Eq. (7). 

Also, for a large center-to-center distance between the 
wires, Sij 3> L, we can expand ai in Eq. ([8|), obtain- 
ing the following expression for the interaction energy 
between two wires 

E,-.iS.,) - , (10) 

This last expression, which we call the dipole-dipole ap- 
proximation, is equivalent to the interaction between two 
dipoles, that is, each wires has been approximated by a 
single dipole. 

In order to investigate the validity of this dipole-dipole 
{E^-^) approximation we calculate the ratio between 



In order to quantify the importance of the interaction 
energy we calculated the ratio between the self-energy 
and the magnetostatic interaction energy between two 
identical nanowires, 

^^ E,r.t{S,,) ^^^^ 
2Eself 

Figure 3 defines the geometry of the two-wire-system for 
which rj = 0.2, 0.1, 0.01 and 0.001. From this figure we 
observe a strong dependence of the interaction energy 
on the geometry of the array. As an illustration, when 
we consider two nanowires with L = 1 /im, i? = 20 nm 
and L/R — 50, if we look for an almost non interacting 
regime, rj — 0.01, 2R/Sij = 0.068, and then the two wires 
have to be at least 590 nm apart. For this geometry 
the interaction energy is about 1% of the self energy. 
However, for the same L and R, if the wires are 58 nm 
apart {2R/Sij = 0.69), the interaction energy is about 
20% of the self energy (ry = 0.2). 
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FIG. 3: 77 values as a function of the geometry of the two- 
wire-system. 



2. Wire array 

We are now in position to investigate the effect of the 
interwire magnetostatic coupling in a square array. Cal- 
culations for the total interaction energy E^^^ (N) of the 
N — n X n square array are shown in Appendix B, and 
lead us to write 

ji-i 

Etrrayin) - 2n ^ (n - (± If (pd) 



(12) 

where + (— ) refers to parallel (antiparallel) magnetic or- 
dering of the nanowires in an array with nearest-neighbor 
distance d, and Eint is the interaction energy between 
two wires given by Eq. ([S]). Note that in an array 5^- is 
a function of d. In the antiparallel configuration the 
magnetization of nearest-neighbor nanowires points in 
opposite directions. Figure 4 illustrates the behavior of 

^arraj/(^) — -^arraj/ ("■) /""^ ^ function of 77, in a fer- 
romagnetic (a) and an antiferromagnetic (b) array. We 
consider an array of identical wires with R — 20 nm and 
L = 1 /im and two different nearest-neighbor-distance d. 
We can see that a large number of wires {N « 10^), cor- 
responding to a sample of w 0.01 mm^, is required for 
reaching convergence of W^„^y{n). However, in view of 
cancellations originated in the different signs of the par- 
allel and antiparallel interactions, the antiparallel config- 
uration converges faster, requiring only the order of 10^ 



wires and a sample of « 1 /im^. 
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FIG. 4: Warrayin) lu a square array of identical wires (R = 20 
nm and L = 1 /xm). We consider different number of wires 
in the array and two nearest-neighbor-distances d. 

We also investigate, for the same wires (i? — 20 nm 
and L —\ /^m), the variation of the asymptotic value of 
^array i''^) ^ functlon of the nearcst-neighbor-distance 
d in a ferromagnetic, (Fig. 5a), and antiferromagnetic 
array, (Fig. 5b). Figure 5 illustrates our results showing 
that in the ferromagnetic array, interaction effects decay 
in an exponential way and extend over long distances as 
compared with d. Figures 4 and 5 agree with conclusions 
from experiments by Nielsch et aUTf who assume that, 
due to the high aspect ratio of the magnetic nanowires 
in an hexagonal array, the stray field interaction extends 
over several nearest-neighbor-distances. 



III. CONCLUDING REMARKS 

We have obtained analytical expressions for the magne- 
tostatic interactions between wires. By expanding these 
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V. APPENDIX A 

In this appendix we calculate the interaction en- 
ergy of two identical nanowires. We start by replacing 
the functional form Mi(r) = MoaiZ in Eq. ([2]), leading to 



Eint ihj) = Mo / A'Uz^^dv 
Jv <Jz 

The magnetostatic potential is given by 



(13) 



, , 1 /■ V • M(r') ^ , 1 /■ n' • M(r') ^ , 

An r - r' At: r - r' 



where V and S represent the volume and surface of the 
wire, respectively. 

Due to the functional form of Mi(r), the volumetric 
contribution to the potential is zero. For the calculation 
of the surface contribution to the magnetostatic poten- 
tial, we use the following expansion on the potential psj 



FIG. 5: Asymptotic values of the interaction energy in a 
square array obtained with A'^ — 10*' (ferromagnetic array) 
and A'' — 10^ (antiferromagnetic array). 
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r - r' 



J2 e""W-^') / dkJ^{kn)J^{kr')e~''^'>-'<^ 



Then, the surface contribution to the potential reads Q 



expressions we investigate first and second order approx- 
imations to the interaction energies showing the range of 
validity of these expansions. When the wires are apart 
distances much more larger that their diameter, the first 
order approximation is valid. We also conclude that the 
dipole-dipole interaction is valid only when 2R/Sij ~ 0. 
The energy expressions lead us to investigate the exten- 
sion of the interwire magnetostatic interactions in the 
array. The number of wires required to reach the con- 
vergence in the energy depends on the relative magnetic 
ordering between nearest-neighbor wires. In the parallel 
array, a big number of wires is required until convergence. 
This result implies that the size of the array is an impor- 
tant factor to be consider when different measurements 
have to be compared. 



oo 2- R 

Uj = Moaj / rf^'e""''^^'^'^ / r'dr' 



m— — oo , 



dkJ.m{kri)Jm{kr')[e 



-k{L-z) _ g-kz-\ 



After some manipulations, it is straightforward to obtain 
the potential at a distance from the axis of the wire 
j-th, Uj 



OO 

RMoaj f dk 



2 J k 





MkR)Jo{kn)[e 



-k(L-z) _ g-fezi 



Now we need to calculate the dipolar field due to the wire 
j-th experienced by the wire i-th a distance Sij apart. In 
this case 
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n = + si - 2rSij cos(</. + /3) , 

where /3 is an arbitrary angle and r defines a particular 
point in wire i. Then, we include in the potential the 
following expansion [25j 



[k^r^ + SI - 2rS,j cos((/) + /?)) 

oo 



which lead us to obtain 

oo 

RM^dj f dk 
[/,- = ' 



yJi(fci?)[e 



oo 
m— — oo 

By replacing this last expression in Eq. ()13p we obtain 

oo 

|jj2(fci?)Jo(fc^,,)(l-e"'=^) . 



Finally, and changing the variable g = kL, we obtain the 
reduced expression presented in the Eq. ([5]), 



Eint{i,j) 



oo 



which can be used to write the interaction energy in a 
compact form; that is 



n n n n 



1=1 j=l k=l 1=1 

We can reduce the number of summations using the fol- 
lowing rule 



^^5(^-*)-"5(0) + 5](«-p)[.9b) + 5(-p)], (16) 

2=1 fc = l p=l 

which lead us to write 



=1 k=l 



n-l 

+ 2Y,{n-p)f^{pJ-j) 



Then, the interaction energy, Eq. ([TS]), reduces to 



J = l /=1 

n — 1 n n 

p=i j=i /=i 

Using again the rule (|16p . we can reduce the double-sums 
in Eq. (I17p . obtaining 



VI. APPENDIX B 

The total interaction energy of the N — n x n square 
array can be written as 



n n n n 

^EEEE(±i)^'=- 

:1 fc=l /=1 



i=l i=l fc=l /=1 



where -I- (— ) refers to parallel (antiparallel) nearest- 
neighbor magnetic orientation of the nanowires in the 
array. Here Eint (0) = 0, avoiding the self-interaction of 
the wires. For simplicity we define the following function 



EU^)='^f^{Q,Q)+nY,{n^q)f^{Q,q) 



n-l 



9=1 
Ti — 1 n — 1 



+n J2 in-p)f^{p, 0) + 2 J2 Y.^n~p){n~q)f^{p, q) . 

p=l p=l g=l 

From (|14p we know that /*(0, 0) = 0, which lead us to 
finally obtain 



Earrayi""^) 



2nJ2in-p){±irE,ntipd) 
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2 E E(" -?')(^^ - g)(±l)P-'^.„t (d^f^^^ 



--1 q=l 



f^{p,q) - (±l)f-'^,„t(dv^^+^) , (14) 
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